Colossal spin fluctuations in a molecular quantum dot magnet 
with ferromagnetic electrodes 
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We study electronic transport through a magnetic molecule with an intrinsic spin S coupled to 
«j ■ two magnetic electrodes, in the incoherent regime. The molecule is modeled as a single resonant level 

' with large Coulomb repulsion (no double occupancy). The molecular spin is isotropic and it interacts 

, with the electronic spin through an exchange interaction. Using an alternative method to the usual 

master equation approach, we are able to obtain analytical formulas for various physical quantities 
of interest, such as the mean current and the current fluctuations, but also the mean value of J z -the 
z component of the total spin on the molecule- and its fluctuations. This allows us to understand 

■ how the electronic current between the magnetized electrodes can control the polarization of the 
molecular spin. We observe in particular that the fluctuations of J z reach unexpectedly high values. 

(N : 

i r I. INTRODUCTION 

13 : 

Molecular spintronics is at the convergence of two recent and rapidly developping fields. On the one hand, molecular 
I . electronics^ where individual molecules are connected to electrodes of different nature, and the effect of the various 
molecular degrees of freedom on the electronic transport can be studied and possibly engineered^ On the other hand, 
spintronics, where the focus is placed on the electronic spin as a new useful degree of freedom. Magnetic molecules - 
molecules having an intrinsic spin, possibly large 3 - - play, of course, an important role in molecular spintronics. Several 
transport experiments have been performed in the past years on such magnetic molecules - specially on molecular 
magnets like Mn%% derivatives, which are molecules with a large spin anisotropy, tending to align the spin along 
an easy-axis^ Theoretical calculations on transport in the incoherent regime for these molecules have been done, 

■ specially in the case where the electrodes have magnetic properties^ Other magnetic molecules are spin isotropic, 
and some work has been done already to study electronic transport with such spin-isotropic magnetic molecule^ In 
Ref. Q, the full counting statistics (FCS) for such a molecule placed between non-magnetic electrodes has been 

i obtained. However, to the best of our knowledge, no work has been devoted to the electronic transport between 

ferromagnetic electrodes through a spin-isotropic magnetic molecule; one of the aims of this paper to study this 
i problem. Experimentaly, transport through such a spin-isotropic magnetic molecule can be obtained for example 
with a magnetic atom trapped inside a Ceo molecule, which is placed between two electrodes (see Fig. [T])j£ 

This paper focuses on the transport through a spin-isotropic magnetic molecule, in the regime of weak coupling to 
the leads. The primary goal is to compute the current, its zero-frequency fluctuations and more importantly to analyze 
the fluctuations of the total spin on the molecule. This is an important issue because of the mutual influence between 
the electronic current passing through the molecule and the molecular spin. We consider ferromagnetic electrodes 
with collinear alignments of electrodes (parallel or anti-parallel, or situation with only one polarized electrode).— The 
system displays a rich variety of behaviors: on the one hand the suppression of the current by spin-blockade and, on 
(*~*) the other hand, unusually large fluctuations of the molecular spin. 

We are considering the temperature regime in which successive tunneling events through the dot are all incoherent 
(incoherent tunneling regime), and describe them as a Markovian process. Such a situation is realized at a temperature 
O which is much higher than the typical energy scale determined by T, i.e., HT <C fcg6. We model the molecule as 
a quantum dot with a single resonant level, with infinite Coulomb repulsion (no double occupancy). The molecular 
spin S , and the electronic spin on the dot level a, interact through an exchange interaction, —J ex S ■ a. The dot level 
is, therefore, split into two levels, corresponding to two eigenvalues of the total angular momentum, J = S ± 1/2. 
For simplicity, we will work at temperatures much smaller than the applied bias voltage, at which electron transport 
happens only in one direction, and the bias window is infinitely sharpJ^ We also choose the chemical potentials of 
the electrodes such that only the J = S + 1/2 spin sector is in the bias window, and only states in this energy level 
take part in the transport (see Fig. [1]) . 

The basic mechanism at work here is the exchange of spin between the itinerant electrons and the molecule: an 
incident spin- up electron, for example, can be collected as a spin-down electron (spin- flip), if the molecular spin has 
its polarization along the reference axis increased by one. Magnetic electrodes, with different densities of states for 
spin up/down electrons, leading to spin-depedent tunneling rates: R = Tl,r(1 +pl,r)/2, R = Tl,r(1 —pl,r)/2, 
can thus induce polarization of the molecular spin (pl,r, £ [—1, 1] is the polarization of the L/R electrode, and Tl r 
is the tunneling rate to the L/R electrode). 
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FIG. 1: Energy diagram of the system: because of the exchange interaction between the molecular spin S and the electron 
spin, the occupied dot level is split between the J = S — \ levels (with J z ranging from — (S — j) to [S — §)) and the J — S + ^ 
levels (with J z ranging from — S — | to S + |), with a spliting ~ J ex S. We chose the chemical potential of the electrodes such 
that the J = S + \ levels only are in the bias window. The inset shows a schematic view of a possible experimental realization: 
a magnetic atom with spin S trapped inside a Cm molecule, which is placed between two electrodes. 



The standard method in the incoherent tunneling regime is to use master equations. This approach has been suc- 
cessfully implemented to compute numerically the current and the noise through such molecular systems*^ Recently, 
analytical results were obtained for the full counting statistics in the case of non-magnetic electrodes^ We use here 
an alternative method which allow us to obtain analytical formulas for the case of magnetic electrodes, for the current 
I, the charge Q, the total spin J z and the fluctuations of these quantities. This method has been introduced by Ko- 
rotkov, for computing numerically fluctuations in the single-electron transistor^ It uses a Langevin approach, where 
the transport process is seen as random sequential jumps between neighbouring system states. We have extended this 
method, in order to obtain analytical results for the present problem. 

The paper is organized as follows. In section II we present a concise but self-contained explanation of the method. 
Section III presents and discusses the results we have obtained in the case of a molecular spin S = 1/2, for the mean 
current and its fluctuations, and for the z component of the total spin J z and its fluctuations. Section IV discusses 
how the results are modified in the case of higher spins. Another method of calculation which can be used to obtain 
the same analytical results is shortly explained in section V, and section VI gives the conclusion. A few appendices 
contain some lengthy formulas, and analytical results for molecular spin 5 = 1. 



II. THE SEGMENT PICTURE 



We give a short self-contained derivation of the method, only stressing the points which are different from the 
original work»ii 



A. General formulation 



The time evolution is divided in terms of segments: a segment £ is defined as a scries of random processes which 
begins with a reference state and finishes with the same state. This reference state is arbitrary, and all the physical 
quantities are of course independent of the choice of this state. As the time evolution is given by a Markovian 
series of random transitions, two different segments are totally independent, and any time integral used to compute 
average or fluctuations can be written in terms of average over the segments. In our model, a state \a) can be 
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characterized by the occupation number of the dot level Q, the total angular momentum J and its z-component J 2 , 
i.e., \a) = \Q,J,J Z ). A segment of length M, starting and finishing with state ao is thus defined by the sequence 
ckq > ct\ ^ &2 — * ' ■ ■ — * cx-M-i — * a 0i an d by the duration of each step. The total duration of the segment is 
r[C] = Sm=o Tm ' wnere T m is the time the system stays in the state a m . 

Considering a random variable X(t), we will compute its average X over the measurement time, i.e., 



X 




dtX{t) 



(1) 



and its fluctuations Sxx- 



S 



xx 




dt (X(t) 




(2) 



Here (• • • ) represents a statistical average over random Markovian process, and we take the measurement time T 
sufficiently larger than all other time scales of the system. For each physical quantity X, we define a function X[Q 
which gives the time integral of this quantity over a given segment e.g., for J z , we define, 



M-l 

^ ^ Jz,m ' 
m=0 



(3) 



where J z , m is the value of J z in the state a 



We have similar expressions for all the other physical quantities which 
have a fixed value in a given state a m . For the current operator /, as its time integration gives the transfered charge, 
we need to define the function k[Q which is the number of electrons tranfered from the left to the right electrode 
during segment C- With these functions, the time integral of X{t) in Eq. (|TJ) over the measurement time T can be 
decomposed into contributions from N successive segments, {(i, C2, ■ ■ • , £,n}'- 



rp TV 

/ dtX(t)^Y,X[( n } 
Jo 1 



(4) 



Taking it also into account that the different segments are independent, one can thus replace the statistical average 
with an average over the segments; 



X = 



I 



N 

J 1 ^ / 



<*[&>]> = 



(X) 



(5) 



where (r) = T/N is the mean duration of a segment. A segment C occurs with a probability P[(], giving (X) 
For the fluctuations of X one first notices that 



dt (X(t) — X) 



N N 



£ J2 (*Kn] -XT[C n ]) {X[C n >\ - Xr[C«']) 



(6) 



n— 1 7 



We note that £ n and £ n i are different segments of the same measurement, i.e., the set {£} = {Ci>C2)"'" iCjv} is 
common for the sum over n and n', and the average is an average over the different segments sets {£}. As different 
segments are independent, the terms with n ^ nl can be written as a product of two averages, which clearly vanish, 
i.e., (^[Cn] — Xr[( n ]) — 0. As a result, one finds, 



'xx — TP, 




dt (X{t) 



(X 2 ) + (r 2 ) 





(7) 



We can thus express the fluctuations Sxx in terms of the averages over the segments, such as, (r), (X), (Xt), etc. 
We note that these are averages over either linear (such as, (r), {Jz)), or quadratic (such as, {J 2 )-, {Jzt)) functions. 
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Note that centered moments of higher order (for example the third centered moment ~ 1/T((J Q dt(X(t) — X)) 3 )) 
cannot be easily obtained with this segment technique. Indeed, the equivalent of Eq. ([7]) for a higher moment has 
terms involving averages over the last, incomplete, segment of the time interval [0, T]. The contribution from this 
incomplete segment is negligible (~ 1/T 1 ) for the first and second moment, but is important for higher moments. As 
this last segment is incomplete, its statistics is different from the statistics of standard segments and cannot be easily 
computed. 

As is emphasized in Ref. [ll|, the average over the segments can be done in two steps: 

{X) = (X) lfl = Y,{X[t\)M& (8) 

c 

The first average is carried out over the durations tq,...,tm-\ of an arbitrary sequence of states «o &i —>•••—» 
otM-i — * «o- This first average, which we denoted in Eq. ([8} as (. . .)%, is easy to perform, as r m is given by a 
Poissonian process, with a rate T m , and we leave further details to Appendix A. The second step is an average over 
all the possible sequences of states, with the correct probability P2K] f° r each sequence. As we will show below, it is 
possible in our case to describe the whole set of sequences, and this second average can also be performed analytically. 



B. Construction of all possible segments — case of molecular quantum dot magnet 

In order to perform the second average, we must identify the whole set of possible sequences for the magnetic 
molecule system. As we work in the limit of strong Coulomb blockade, the dot level can be occupied at most by 
1 electron, so we have Q = or Q = 1. For the empty dot, the total spin is simply given by the molecular spin, 
and specifying the z component of the spin determines the state completely (S z £ [-S, S]), so an empty dot state is 
given by |0, S z ). For the occupied dot (Q = 1), the total spin is obtained by the addition of the molecular spin and 
the spin of the electron occupying the dot, which gives J = S ± 1/2. These two values of the total spin correspond 
to two levels of the system, separated by an energy of order J ex (the value of the exchange coupling between the 
spins). As explained before, we decide here to work in the ferromagnetic case, where the lower level is the one with 
,/ = S + 1/2, and with the chemical potentials of the electrodes placed so that only this lower level is in the bias 
window; then the level J = S — 1/2 plays no role in transport (it cannot be populated) and can be forgotten. The 
occupied dot has thus a total spin J = £+1/2, and specifying the z component again determines the state completely 
( J z E [—(S + 1/2), S + 1/2]), so an occupied dot state is given by |1, J z ). 

For the reference state (which is the initial and final state of each segment), we choose an empty dot with spin 
maximally polarized along the z axis: \0,S Z = +S). From this state, there are two basic sequences where a single 
electron is tranfered from the left to the right electrode: 

• (A) |0, S z = +S) —*\l,J„ = S+ 1/2) -» |0, S z = +S) 

• (B) |0, S z = +S) -» |1, J z = S - 1/2} -» |0, S z = +S). 

The two sequences (A) and (B) are the two simplest ones. Clearly, the sequence A cannot be extended further, as 
the spin of the intermediate state is maximal. On the contrary, the sequence (B), which we will call a basic sequence 
B, can be extended by adding a subsequence starting and finishing at the intermediate state \1,J Z = S — 1/2), and 
going only to lower values of S z . One can add 

• ix times the subsequence v\ 

|1, J z = S - 1/2) -» \0,S, = S - 1) -> |1, J Z = S- 1/2) 

• 12 times the subsequence vi 

|0, S z = S - 1) -» |1, J z = S - 3/2) -» |0, S z = S - 1) 

• iis times the subsequence 

|0, S z = -S) -> |1, J z = -S - 1/2) -> |0, S z = -S) 

Combining the 4S 1 subsequences to the basic segment B, with arbitrary repetition of each subsequence, one can 
construct all the possible sequences of type B (see also Fig. [5]). We use the notation B*(ii,i 2 , ■ ■ ■ , i^s) to represent 
the type-B sequence composed of the basic sequence B combined with i\ times the subsequence vi, 12 times the 
subsequence V2, etc. The basic sequence B is then simply B*(0, 0, . . . , 0). The whole set of type-B sequences, plus the 
simple sequence A, give all the possible sequences. 
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FIG. 2: Construction of all the possible sequences - type A and type B+extensions. A sequence starts and ends in the reference 
state, Q = 0, S z — S. Sequence A (solid line) correspond to the tunneling of a spin-up electron from the left to the right lead. 
Sequence B (dashed line) to the tunneling of a spin-down electron from left to right lead. Sequence B can be extended by 
attaching subsequences starting and ending from the state Q = 1, J z = 5—1/2 (thin dashed lines), forming longer sequences 
where several electrons are transmitted, and where the molecular spin goes through intermediate states with S z < S and 
J z <S- 1/2. 



The probability with which a given sequence occurs is given by the product of the probabilities of the transitions 
forming that sequence. The probability of a transition is given by the transition rate divided by the total transition 
rate of the initial state. One must here distinguish between transitions starting from an empty dot (Q — 0) and 
transition starting from an occupied dot {Q = 1). An empty dot (Q = 0) with a molecular spin S z is subject to two 
types of transitions: tunnel-in of an electron from the left electrode with either up or down spin, which brings the 
dot to the state J z = S z + 1/2 or J z = S z — 1/2. The transition rates Fq_ (S z ) and Fq =0 (S z ) for the above two 
processes can be calculated, using Fermi's golden rule and Clebsh-Gordan coefficients: 

r± r<? \ — F 1± plS + i±s z 

Using the tunneling rates ^q =0 (S z ), one can express the probability Pq =0 (S z ) with which the system jumps onto 
either of the two final states: 

Fq =0 (S z 



P ^ {Sz) -rX^sV+rL{s,y 



For an occupied dot (Q — 1), with the molecular spin state J z , tunnel-out of an electron with either spin-down or 
spin-up brings the dot to the state S z = J z + 1/2 or S z = J z — 1/2, with the rates 

r Q=1 (J*)-r fl — 25 + 1 , (ii) 

giving the probability 

^ l(j2) = ^W + r Q ,(J z )' (12) 

Using Eq. (|10ll2p . one can express the probabilities of sequences A and B: 

P2[A] = P+ =q (S z = S)xPq =1 (J z = S +1/2), 

P2[B] = P Q=Q (S Z = S)x P+ =1 (J Z = S- 1/2). (13) 
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The subscript 2 is to recall that these probabilities are associated with the second average (average over different 
sequences). Similarly, extensions starting from an occupied dot or an empty dot occur with the probability, 

PoW = P Q ^(Jz) x P^ (S Z = J z — 1/2) 

P C (S Z ) = P Q=0 (S Z ) x P+ =1 {J Z = S Z - 1/2) (14) 

where the subscript o and e stand for occupied and empty. The probability of the sequence B*(ii, i 2 , ■ ■ ■ , i<is) is then 

P 2 [B*(H,i 2 , us)} = Cl +i2 _, . . . C^_ 1+i4S -i ■■■P% S P2[B], (15) 

where P2/-1 = Po{J z = S — (21 — l)/2), pn = P C {S Z = S — I). In this expression, the combinatorial factors 
{C\ — — count the number of different sequences corresponding to the set (ii,i 2 , ■ • ■ ,*4s) because of the 

possible permutations of the subsequences. There are C| n _ i+i _ 1 different possibilities to "attach" the ii subsequences 
at level I to one of the intermediate states of the subsequences at level I — 1. The second part of the expression 
is simply the product of the probabilities of all the subsequences. 

Together with some explicit formulas for the first average (see Appendix A), Eqs. (|13I15[) allow us to evaluate the 
averages over different segments ( appearing in Eq. ([8|). Evaluating all such averages, one finally finds the explicit 
formulas for X and Sxx, formally written as Eqs. ([5|) and ([7|). Analytic results for such quantities are listed in Tables 
I and II, as a function the polarization p of the electrodes. Some of such examples are also shown in Appendix B, 
along with some intermediate steps in a specific case of S — 1/2. Note that we are able to obtain such analytical 
formulas owing to identities involving summation on the binomial factors, like: 

00 

3=0 v ' 

In practice, the calculations are quite lengthy, but results are easily obtained using a symbolic computation software. 
In Appendix B, it is also shown that our analytic results are consistent with the cumulant generating function, 
obtained in Ref. Q, in the limit of non-magnetic electrodes: p — > 0. 

III. RESULTS 

In this section, we present the results we have obtained for the various quantities. We will focus on the mean 
current / and the current noise Sn, on the mean charge on the dot Q and its fluctuations Sqq, and on the mean 
value of z-component of the spin on the dot, J z and its fluctuations Sj^j z . In order to see the effect of the leads 
magnetization, we will show all these quantities as a function of the polarization of the leads. For simplicity, we choose 
to have a single parameter for the leads polarizations, and we have chosen four representative cases. In the first two 
cases, the two electrodes are magnetic. The absolute value and the direction of the polarizations are the same in the 
two electrodes, but are either parallel (Pl = Pr = p, case noted P), or anti-parallel (Pl = —Pr = p, case AP). In 
the two remaining cases, only one of the electrode is magnetic (with a polarization p), while the other one has no 
magnetic property. The polarized electrode can be either the left one, which is the source electrode (case LP), or the 
right one, which is the drain electrode (case RP). 

Table U gives the anaylitic formulas we have obtained, in the four different cases for the leads polarization, for a 
molecular spin 5 = 1/2 (Appendix [Cl shows a similar table for the case S = 1; we do not show any formula for a higher 
spin S as they become too lengthy). For simplicity, we have chosen equal bare transition rates for the left and right 

electrodes, T ® = = 1 (formulas with general transition rates can be obtained easily with the same method). As 
can be seen on the table, each analytic formula is given by a fraction of two polynomial in p, and containing only 
of even powers of p (except for J z where there is an additional factor p) . The orders of these polynomial, and their 
coefficient, increase when the molecular spin is increased (see Appendix [C]) 

A. Mean current and current fluctuations 

The behavior of the mean current as a function of the leads polarization, for the case of the molecular spin S = 1/2, 
is shown on the left panel of Fig. [3] The simplest case is the case of parallel polarizations (P, dotted curve): the mean 
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TABLE I: The analytic formulas for the case of a molecular spin S = 1/2, when only the triplet state of the occupied dot 
(with a spin 5 = 1) lies in the bias window. The first column shows the type of polarizations in the leads : P for parallel 
(Pl = Pr ~ p), AP for anti-parallel (Pl ~ —Pr = p), LP (Pl = p, Pr — 0) for left lead polarized only and RP (Pr — p, 
Pl = 0) for right lead polarized only. The second column shows the quantity whose analytical formula is given in the third 
column 

current is then constant (value 3/10), and the polarization of the leads has thus no effect. Note that the fact that the 
current is constant is related to the choice we made for the density of states py±j — (1 ±pj)/2 (j = L,R). Indeed, 
we see that the total densitiy of states (spin up + spin down) in each electrode is constant. 

Consider next the anti-parallel case (AP). We see (full curve on Fig. [3]) that the mean current has a maximum for 
zero polarization, and decreases to zero when p approaches ±1. This behavior can be understood simply: when \p\ is 



8 




FIG. 3: The mean current I (left panel) and its flucutations Su (right panel) for the case of a molecular spinS* = 1/2, for the 
four different cases of electrodes polarizations. The inset in the right panel shows the behavior of Su in the parallel case (P) 
on a larger scale. 

large (let us take for example p close to 1), the electrons coming from the left electrode have preferentially a spin up, 
while the electrons going to the right electrode have preferentially a spin down; the transport of such an electron from 
the left to the right electrode implies thus a flip of the electron spin, and thus an increase of 1 of the z component 
of the molecular spin. However, this will lead quickly to a maximally polarized molecular spin, for which such a spin 
exchange will be impossible. The only processes contributing to transport will then involve the electrons with a low 
density of state (spin down in the left lead or spin up in the right one), for which the current goes to when \p\ goes 
to 1. This behavior is a case of spin blockade: for \p\ = 1, an electron is blocked on the dot because its spin does not 
fit the collecting electrode spin. 

The case where the right electrode only is polarized (RP, dash-dotted curve on Fig. [3]) is similar. For p = 1, the 
system reaches a state where the z component of the molecular spin is maximally negative (S z = — S for the empty 
dot). In this state, spin-up electron can tunnel from the left electrode to the right electrode without any splin-flip, 
but as soon as a spin-down electron is tunneling from the left electrode, it is blocked on the dot because it cannot flip 
to a spin-up electron and tunnel to the right electrode. Because of the very large Coulomb repulsion on the molecular 
level, the presence of this spin-down electron forbids any further transport of spin-up electron. This case is thus also 
a case of spin blockade, as in the AP case, but the Coulomb repulsion on the molecular level plays here a central role. 
The decrease of the mean current with p is a bit slower than in the AP case, as the process where an electron tunnel 
without any spin flip, and without involving small densities of states, is always possible when \p\ < 1. 

Finally, the case where the left electrode only is polarized (LP, dashed curve on Fig. [3]) has a totally different 
behavior. There, the current is slowly increasing when \p\ increases, and there is no spin blockade. Indeed, as the 
density of states of both spins are important in the right electrode, it is always possible to have transport of an 
electron without spin flip, and involving large densitites of states. 

Let us now consider the zero-frequency current fluctuations, shown on the right panel of Fig. [3J Again, we observe 
very different behavior for the four different cases of electrode polarizations. At p — 0, the value of Su is 39/125 = 
0.312. In the parallel case (P, dotted curve), Su increases with p, and diverges as (1 — p 2 )^ 1 as \p\ approaches 1 
(see the inset in the figure). On the other hand, in the anti-parallel case (AP, full curve), the current fluctuations are 
nearly maximal at p = (with small shoulders near \p\ = 0.3), and decrease to as \p\ — » 1. This huge difference in 
behavior can be understood using the segment picture; let us take for example p close to 1. In the parallel case, the 
most probable process contributing to transport is simply the transport of one spin-up electron, without any spin flip, 
and thus without exchanging angular momentum with the molecular spin. The most probable segment is thus a very 
short one, with a single electron transfered. However, an exchange of angular momentum (spin flip for the electron, 
and modification of the z component of the molecular spin) can happen with a small probability; when this happen, 
the system will then transfer again a very large number of electrons without any spin flip, and it will take a very 
long time before the z component of the molecular spin recover its inital value. There is thus a small probability to 
have a very long segment, with many electrons transfered - the smaller the probability, the longer the segment. This 
presence of rare but arbitrary long segment when p goes to 1, in a "background" of very short segments, explains the 
divergence of Su in the parallel case. The situation is different for the anti-parallel case: for p close to 1, the molecular 
spin is with a high probability in a maximaly polarized state (S z = S for the empty dot). The most probable process 
is again the transfer of a single electron without spin flip (this produces a low current as it involves a small density of 
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FIG. 4: The mean value of the z component of the molecular spin, J z (left panel) and its flucutations Sj z j s (right panel) for 
the case of a molecular spin S = 1/2, for the four different cases of electrodes polarizations. 

state in one of the electrodes). There is again a small probability of a spin- flip, which will bring the molecular spin 
in the state S z = S — 1. However, at this point, the most probable process (involving large densities of state in both 
electrodes) tends to bring the molecular spin back to the S z — S state. The resulting segment is thus also short, with 
2 (or at most a few) electrons tranferred. As this probability for such longer segments goes to zero when p — ► 1, and 
as the length of these segments is quite short, we understand why the current fluctuations Su go to when \p\ — > 1. 

In the case where the right electrode only is polarized (RP, dash-dotted curve), the current fluctuations also go to 
as \p\ — > 1. Note however the presence of broad shoulders, with a maximum of the fluctuations near \p\ — 0.6. Finally, 
in the case where the left electrode only is polarized (LP, dashed curve), the fluctuations have a behavior similar to 
the one of the mean current, with a slow increase when \p\ increases. 

B. J z and the J 2 fluctuations 

In many works about electronic transport in the incoherent regime, the emphasis is put on the statistics of the 
electronic current, and little attention is given to the statistics of other quantities (see however Ref. [l|). Here we 
study the statistics of the total spin of the molecule (i.e. its z component J z ), which gives us precious information on 
the impact of the electronic current on the molecular spin. 

The results for the mean value of the z component of the molecular spin, J z , and the J z fluctuations Sj z j z> are 
shown on Fig. 2] Note that these two quantities involve both the molecular spin when the dot is full (J z ) and the 
molecular spin for the empty dot (which is noted S z , but as the intrinsic molecular spin S is also the total spin for 
an empty dot, J z reduces to S z for an empty dot). The mean value J z is an important quantity, as it shows how 
the current through the molecule is changing the polarization of its spin (as without any current, one has simply 
J z — 0). The fluctuations Sj z j z show how the molecular spin fluctuates around its mean value, it gives thus precious 
information on how precisely one could control the molecular spin polarization by applying a current. 

In the parallel case (P, dotted curve on the left panel of Fig. QJ, the mean value J z is simply for all p. In the 
anti-parallel case (AP, full curve), J z is an non-linear odd function of p, going from to 1 for p going from to 1. 
This behavior is easily understood. For p > for example, the spin-up electrons have a larger density of state than 
spin-down electrons in the left electrode, and it is the opposite in the right electrode. The system thus favors the 
transport of a spin-up electron from the left electrode into a spin-down electron in the right electrode (compared to 
the process with the spins exchanged), and this process increase J z by 1. On average J z will thus be positive for 
p > 1. For p — 1, we have seen that current is because of spin blockade. In this case, the system is frozen in the 
state where the dot is full, with J z = 1, hence one has J z = 1. 

The behavior is quite similar in the case where the left electrode only is polarized (LP, dashed curve): J z is an odd 
function of p, positive for p > 1. There are two main differences with the anti-parallel case. First, the slope at p = 
is smaller (it is 1 for LP, and 2 for AP). Next, J z does not reach 1 (but 5/6) for p = 1. This is because the current is 
non-zero even for p = 1 (no spin-blockade), and thus the molecular spin oscillates from the value J z = 1 (filled dot) 
and S z — 1/2 (empty dot). Finally, the case where the right electrode only is polarized (RP, dahs-dotted curve) has 
an opposite sign, and the behaviour is simply linear, with J z = —p. The fact that the sign is the opposite from the 
sign of the other cases ( J z — —1 for p — 1) is due to the fact the tunneling of a spin-up electron to the right electrode 
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(dominant for p > 1) is done either without spin flip (if the electron coming from the left electrode is also a spin- up 
electron) , or with a spin- flip which decreases J z by 1 (if the electron coming from the left electrode is a spin-down 
electron). The fact that J z reaches the value — 1 for p — 1 is again due to spin blockade, as in the anti-parallel case. 

Let us now consider the fluctuations of J z around its mean value, Sj z j z (right panel of Fig. 2]). Note first that the 
value of Sj z j z for p = is 88/5 = 17.6. This value, which determines the overall scale of the fluctuations, is extremely 
large, and is discussed in more details below. In the parallel case (P), Sj z j z is increasing as \p\ increases, and has the 
same (1 — p 2 )~ l divergence as the current noise (for the same reasons). In the other three cases, Sj z j z is maximum at 
p = 0, and decreases as \p\ increases, with a much broader shapes for the cases with only one electrode polarized (LP 
and RP) compared to the anti-parallel case (AP). In the anti-parallel case, and in the case where the right electrode 
only is polarized (RP), Sj z j z is for \p\ — 1, because of spin-blockade (current is zero). But in the case where the left 
electrode only is polarized (LP), Sj z j z = 4/27 for \p\ = 1 because there is no-spin blockade: for p — 1 the molecular 
spin oscillates between the values S z — 1/2 and J z = 1, leading to these non-zero fluctuations. Finally, we note that 
the behavior of Sj z j z is very similar in the two cases where there is only one electrode polarized (LP and RP); this 
is quite remarkable, as the currents (and the current noises) in these two cases have a completely different behavior 
(see Fig. El). 

As said above, the scale of these J z fluctuations is very large, with a value Sj z j z = 88/5 for p — 0. A natural 
normalization of these fluctuations, to take into account the value of the molecular spin, is to divide by (25 + l) 2 
(note that h = 1). Here, 5 = 1/2, which gives a normalized value of 22/5. This value has the dimension of a time, 
and should be compared with a typical time in the system. Here, the natural time is just the inverse transtion rate 
1/r^^j = 1, which gives the scale of the time to transfer an electron. The value of the J z fluctuations at p = is 
quite larger than this time scale. Comparison with the value for the charge fluctuations Sqq, which is 96/125 ~ 0.77 
(see table [J) shows also that the normalized J z fluctuations for p = are extremey large. Even if these fluctuations 
decrease with increasing \p\ (except in the parallel case), they remain quite large when \p\ is not close to 1. One can 
thus speak of colossal spin fluctuations, and this implies that it is difficult to control the molecular spin with the 
current, except with electrodes having polarizations p close to 1. 

In this respect, the anti-parallel case is much more favourable than the case where only one electrode is polarized. 
One could for example think to use a setup with only one polarized electrode, to flip the molecular spin by reversing 
the current in the setup. Indeed, reversing the voltage bias will make the system go from the LP case to the RP case. 
If p = 0.5 for the polarized electrode, we see on the left panel of Fig. [4] that J z would change from approximatively 
0.5 (case LP) to —0.5 (case RP) when the bias voltage is reversed. However, as the J z fluctuations are very large (the 
normalized value is ~ 5), it is difficult to say that the molecular spin is controlled. Performing the same with two 
polarized electrodes in the anti-parallel polarization configuration would be more effective: reversing the bias voltage 
is then equivalent to the change p — > — p, and with p = 0.5 it would change ( J z ) from approximatively 0.75 to —0.75. 
The normalized value of the J z fluctuations is then approximatively 1.3, which is much lower than in the previous 
case. 



IV. BEHAVIOR FOR LARGER MOLECULAR SPIN 

In the previous section, we have shown the results obtained in the case of a molecular spin 5 = 1/2. The method 
we have presented is of course not limited to this value of the spin, and we discuss in this section how the results are 
changed when the molecular spin is larger than 1/2. Note that the analytical results for 5 = 1 are given in table [H] 
in appendix [C] The plots one obtains with these results are qualitatively similar to the ones for 5 = 1/2. 

In order to discuss the behavior at larger 5, one should distinguish between the mean current (and the current 
fluctuations) on the one hand, and the mean value of J z (and the J z fluctuations) on the other hand. For the mean 
current, and the current fluctuations, there is very little change as one increases the spin, and the physical explanations 
we have given for S = 1/2 apply for arbitrary spin. This is illustrated on Fig.[5j which shows the mean current I (left 
panel) and the current fluctuations Sn (right panel) as a function of the polarization p, for the anti-parallel (AP) 
configuration, and for the values of the spin 5 = 1/2, 5 = 1 and 5 = 3/2. One can see that both the mean current 
and the current fluctuations decrease a little bit as the spin is increased, with very little change in the p dependence. 
For larger values of 5, the curves will slowly converge towards a "classical curve" , obtained by considering a classical 
(fixed) spin. 

The situation for the mean of J z and its fluctuations is slightly different. First, as J = 5 + 1/2, it is natural to 
normalize the results for the different 5 to compare them; we normalize J z by 5 + 1/2, and the J z fluctuations by 
(25 + l) 2 . The normalized mean of J z is shown on the left panel of Fig [HI for the case of anti-parallel polarizations 
(AP), and for the values of the spin 5 = 1/2, 5=1 and 5 = 3/2. When p — » 1, the molecular spin is maximaly 
polarized, and thus each curve reaches the value 1 for p = 1. However, the slope at p = increases when 5 increases. 
By inspection of the formulas, we see that the slope at p = is given by 4/3 * (5 + 1). The molecular spin is thus 
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FIG. 5: The mean current (7) (left panel) and its fluctuations Sn (right panel), in the anti-parallel configuration of the 
electrodes, for the values of the spin S = 1/2, S = 1 and S = 3/2 as a function of the electrodes polarizations p. When 
increasing the spin S, the curves converge towards a "classical" curve (dotted line) corresponding to a fixed spin. 
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FIG. 6: The normalized mean value of J z , (J Z )/(S + 1/2) (left panel) and its normalized fluctuations S.j z j z /(2S + l) 2 (right 
panel), in the anti-parallel configuration of the electrodes, for the values of the spin S = 1/2, S — 1 and S = 3/2 as a function 
of the electrodes polarizations p. 

more easily polarized when S increases. When the molecular spin becomes large, we expect that it becomes more 
sensitive to the electrodes magnetizations, reaching even for small \p\ highly polarized states. 

The J z fluctuations, Sj z j z1 normalized by (2S + l) 2 , are shown on the right panel of Fig [SJ Two important 
characteristics appear on this figure. First, the maximum value of the normalized fluctuations, for p = 0, increase 
strongly when S increase. This maximum is already very large for S — 1/2 (see the discussion in the previous 
section), but it is still much larger for larger S. At p = 0, Sj z j x = 88/5 for S — 1/2, Sj g j z = 552/7 for 5 = 1, and 
Sj,j x = 2080/9 for S = 3/2, • • • . Secondly, the width of the curves as a function of p decreases as S increases: the 
full width at half maximum is ~ 0.6 for S = 1/2, ~ 0.45 for S=l and ~ 0.35 for S = 3/2. This means that for \p\ 
large enough, the normalized fluctuations decrease when S is increased. For example, on the figure, we see that for 
\p\ > 0.4, the normalized J z fluctuations for S = 3/2 are smaller than those for S = 1/2. We expect this tendency 
to continue when S is increased, with the normalized J z fluctuations for a large spin S having the shape of a narrow 
peak with a very large maximum value. As soon as the electrodes have some magnetization, a larger spin is thus 
relatively easier to control than a small spin S, as the normalized J z fluctuations can be much smaller. 



ALTERNATIVE METHOD OF CALCULATION 



We discuss here shorlty another method of calculation, which can also be used to obtain analytical formulas for 
the averages and fluctuations of different quantities, and which can even give access to analytic expression for the 
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higher moments. It is derived from the master equation approach to the full counting statistics, which was introduced 
in Ref. [13j|. and adapted for a molecular quantum dot magnet (placed between normal electrodes) in Ref.[7]. A full 
explanation of the method can be found in these two references. 

In the master equation approach, a n x n matrix L determines the time evolution of the populations of the n 
different states of the system : 

where pit) is the vector containing the populations. Any off-diagonal element of the matrix L gives the transition 
probability from state i to state j. To obtain the full counting statistics of a given physical quantity, a counting field 
£ is introduced in the matrix L, by making the appropriate replacements of diagonal and non-diagonal elements. The 
eigenvalues of the matrix L(£) then give access to the full counting statistics, as the cumulant generating function 
is simply proportionnal to the eigenvalue A(£) which satisfies lim^o -M£) = 0- The full counting statistics is thus 
obtained by solving the equation: 

A" + U-iiO^' 1 + ■■■ + + /o(0 = 0, (18) 

where the functions /i(£) (i = 0, . . . , n — 1) depend on the matrix L modified by the counting field £. The solution 
A(£) which satisfies lini£„>o A(£) = then gives access to the cumulants C n (X) of the quantity X associated with the 
counting field £: 

Note that the first two cumulants (Ci and C2) are simply the average and the fluctuations which have been calculated 
in the previous sections. In general, the cumulant of order n can be expressed as a combination of the centered 
moments of order < nJ^ T is the measuring time, which must be larger than all typical times in the system. It is in 
general impossible to solve Eq. (fT8|) analytically (except in special cases which can be reduced to small n, as in Ref.0]). 
However, if one is interested in the cumulants up to a finite order n ma x only, then one can expand the function A(£) 
in power of £, and keep only the terms up to the order n ma;r . It is then possible to solve Eq. (fT8|) by expanding all the 
terms in powers of £, and by solving order by order, starting from order 1, up to order n max . Specifically, to compute 
the average and fluctuations (of the chosen quantity) only, it is enough to write A(£) = Ci£ + (C2/2)£ 2 , to develop 
/o(0; /i(0 an d /2(C) U P to order 2 in £, and then to solve Eq. (|18[) first for C\ (terms in £) and then for C2 (terms 
in 6). 

There are two kinds of observable with different types of countings fields. First, charge- like operators, which have 
a given value for each state of the dot (for example the charge Q, or J z ). In this case, the counting field is simply 
introduced by adding +c Q £ to each diagonal element L aa , where c Q is the value of the observable (for example J z ) in 
state a. Secondly, current-like operators, which are associated with transitions between different states (the charge 
current being the main example). In this case, the counting field is introduced by multiplying the off-diagonal elements 
L a p which are associated to transitions contributing to the current by e** . For this second kind of observable, there 
is an additional factor i n on the right hand side of Eq. (fP5)) . 

With this alternative method, we have computed all the quantities shown in the previous sections of this article, and 
we verified that we could indeed recover the same formulas. This method allow to compute quite easily cumulants of 
order higher than 2. We do not provide here a complete exploration of the higher cumulants of the physical quantities 
we are interested in, but as an example FigjT] shows the third and fourth cumulant of the current for the four possible 
cases of electrode magnetizations. We see on these plots that some of the features present for the current fluctuations 
(Figj3]) are more pronounced on the third and fourth cumulant. In the right-lead polarized case (RP), peaks for 
large p are becoming more pronounced, while in the left lead polarized case (LP), the variations of the cumulant as a 
function of p are extremely small. 

VI. CONCLUSION 

In this paper, we have studied the incoherent transport between ferromagnetic electrodes through a magnetic 
molecule with an isotropic spin. The molecule is modeled as a single resonant level, with large electronic interaction 
forbidding double occupancy of the level. There is an exchange coupling between the molecular spin and the electronic 
spin on the molecular level. By extending an original method (which was first introduced by Korotkov in the study 
of noise of the singe-electron transistor—), we have shown that it is possible to obtain analytical formulas for the 
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FIG. 7: The third cumulant of the the current 63(1) (left), and the fourth cumulant of the current C&{I) (right) as a function 
of the electrodes polarizations p, for the four different cases of electrodes polarizations. The inset on the right panel shows the 
curves on a larger scale. 



average and fluctuations of all physical quantities of interest. The idea of the method is to separate the transport 
process in statistically independent segments, and to compute the average and fluctutations using the properties of 
individual segments and then averaging using the statistical distribution of all segments. We have focused on the 
current through the molecule and its fluctuations, and on the total spin of the molecule and its fluctuations. We 
limited the calculations to the case of a temperature much smaller than the bias voltage, and with only the level 
J = S + 1/2 in the bias window, but using the same method it would be quite easy to generalize the results to higher 
temperatures, or to the case where the two levels J = S ± 1/2 are inside the bias window. 

By considering several configurations of the magnetization of the electrodes (both electrodes polarized with parallel 
or anti-parallel polarizations, or only one electrode polarized), we have shown that there is a rich variety of behaviors. 
Depending on the electrode polarizations, the current can decrease or increase with the electrodes polarization; for 
maximally polarized electrodes, it can be blocked due to spin blockade. The current fluctuations can show a non- 
monotonic behavior as polarization is increased. 

With the study of J z (the z component of the total spin of the molecule), we have been able to characterize how 
the electronic transport affects the molecular spin. If the results for the average of J z show, as expected, that it 
is possible to polarize the molecular spin by using magnetic electrodes (non-zero average of J z whose sign depends 
on the sign of the current), the results for the fluctuations of J z show that these fluctuations are very large. This 
shows that it is effectively difficult to control the molecular spin with the current, except with electrodes having 
polarizations close enough to 1. We have shown how these results evolve when increasing the bare molecular spin: the 
polarizability of the molecular spin increases near p = 0, and the fluctuations of J z are more peaked around p = 0. 
These large fluctuations of the molecular spin are of course a direct consequence of the isotropy of the molecular spin. 
The fluctuations would be severly reduced in molecular magnets, where a strong spin-anisotropy is present. 

We have also shown that we can obtain the same results using a different calculation, based on an extension of the 
method introduced by Bagrets and Nazarov 13 to compute the full counting statitics in Coulomb blockade systems. 
In contrast to the segment method, this second method offers the possibility to compute the higher cumulants of the 
physical quantities, and we have shown as an example the results for the third and fourth cumulant of the electronic 
current. This second method gives however less information for the physical interpretation of the results, compared 
to the segment method. 
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APPENDIX A: EXPLICIT FORMULAS I: FIRST AVERAGE — AN AVERAGE IN A GIVEN SEGMENT 

In the segment picture, the statistical average X and the fluctuations Sxx of a given random variable X(t), reduces 
to an average of quantities in a single segment, such as, (X), (t), (AY), ■ • ■ , and its lowest order examples are, Eqs. 
(|5l7p . In these formulas, the average in a single segment is done in two steps, i.e., first a Poissonian average over a given 
segments is taken, and then one takes an average over different segments with suitable weights. One may denote the 
first average for a given segment £ as (• • • )i- The probability with which in a given segment C — (c^ij o^, ' ' ' > a M-i)i 
the m-th state happens to survive during a period r m is proportional to exp[— r m r m ], where l/r m is the average 
lifetime of the m-th state. The probability with which Q is characterized by a set of transition times (to,ti, • • • , tm-x) 
thus obeys to a Poissonian distribution, 

Pi(t ,ti, ■ ■ ■ ,tm-i) = r r! ■ ■■r M ~ie- r ° Ta e~ r ^ ■ .. e - r ™-^>-\ (Al) 

and in that case, one can rewrite the probability P[(] of the segment £ as P[(] — Pi(tq,ti, ■ ■ ■ , tm-i)-P2[C]- Then, 
together with the explicit definition of the first average, 

(*[C]>1 = / dT dn ■ ■ ■ dTM-lX[£]Pl(T , 71, • • • , 7M-l), (A2) 

one indeed arrives at Eq. p. In Eq. (|A"2|) . X[(] is, e.g., J Z [Q defined as Eq. j3]). As for the second average, Section 
II B demonstrates how to construct explicitly P2[C] i n the case of molecular quantum dot magnet. 
At linear order, the role of (• • • )i is nothing more than a replacement r m — > (r m )i = 1/T m , e.g., 

M-l M-l Jz 

mc])i=Y. Jz m^)i=H^- (A3) 

r, n m 

At second order, one may also use the following mathematical trick^ Let us attempt to calculate, e.g., {j7z[C] r [C])i : 

I M-l M-l \ I M-l \ 

{Jz[(\r[Qx = I ]T J ^ra E T ™') = ( E JZ m^ra) 2 + E J^T m T m ,j 
\ ra=Q m'=0 I 1 \m=0 m=£m' I ^ 

M-l 

= zl J rn.{( T in) 2 )l + ^ J m( T ™.T m >}l 
m=0 m^m' 
M-l 

= V j z — + V j z — — 

/ 4 in p2 / j m -p p 

n rn / , L m L m' 

m—0 rn^m' 

M-l 

= E + a[C])iMC])i (A4) 

At the second identity, we divided the double summation into diagonal and off-diagonal parts. To give some concrete 
examples, in the case of segments ( — A and £ = B*(ii, i 2 , ■ • ■ , lis) defined in Section II, the final expression reads 
explicitly as, 

WMB-l), = + J' 1/2 p (A5) 



2s 

E*2i-1 



S -1 + 1/2 S-l 




i=i \ I 1 j z =s-i+i/2i L L s x =s-il 
+ (J z [B*])i{t[B*])i, (A6) 

where we used an abbreviated notation, Ig~° = Tq = o(S z ) = Fq—q(S z ) + Fq=q(Sz), and Ij _ 1 = Tq = i(J z ) 
t q=i( j z) + r Q=i( J z)- Substituting Eqs. (|A5IA6|) into Eq. © in Section II A, and together with P 2 [A] and P 2 [B*} 
constructed in Section II B (Eqs. (|13ll5p ). one finds (J z t) appearing in Eq. for X = J z . Evaluating other 
averages in a single segment, such as (r), (J z ), (t 2 ), and {J z ), one finally finds the expression for Sj g j, = 2/j, 2 [J z ]. 
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APPENDIX B: EXPLICIT FORMULAS II: SECOND AVERAGE — AVERAGE OVER DIFFERENT 
SEGMENTS AND CONSISTENCY WITH THE FCS GENERATING FUNCTION 



To illustrate how to evaluate the average over different segments, let us give here some explicit formulas, which 
typically appear in the calculation. For simplicity, we consider here only the case of S = 1/2, and parallel (P) or 
anti-parallel (AP) spin alignment of the electrodes: Pl = Pr = p for P, and Pl — —Pr — p for AP. Using Eq. ([8]) 
in Section II A, and some explicit formulas in Appendix A and in Section II B (expressions for -p2^4] and P2[B*], in 
particular, i.e., Eqs. (|13ll5p ). one finds the average duration of a segment as, 

2 3r L + 2T R 
3+p T l Tr 

, v 2(3 + io P 2 + 3p 4 )r L + 4(i-p 4 )r fl 

{t)ap = (i +P )3( 3 -2 P -p»)r L r fl ■ (B1) 

Here, we used a slightly different convention from the body of the paper, so that we can compare our results directly 
with that of Ref. 0. The convention here is T\ R = Tl{1+pl,r) and T l L R = T L (1 -pl.r)- Note that Eqs. (|Bl~j) are 
not symmetric functions of p, and this reflects the choice of the reference state. The average net charge on the dot is, 

<0\ 6 (O) 2(3 + 10p 2 + 3 P 4 ) 

(Q>p -3T? {Q)Ap - (i+p)3(3-2p- P »)r Ji ' (B2) 

no longer a symmetric function of p. On the other hand, the average charge normalized by the average duration of a 
segment, which is the physically measurable average charge, is a symmetric function of p, indepedent of the choice of 
the reference state: 



(t) p ST l + 2T 



r 



n _ (Q)ap = (3 + 10p 2 + 3p 4 )r L 

^ AP (t)ap (3 + 10p 2 + 3p 4 )r L + 2(l-p 4 )iV 1 6) 



As for the current, one has to evaluate similarly, (fc)p or (k) ap to find, 

(k) P 3r L r R 



Ip 



(t) p sr L + 2T R ' 



j (3ap (1-p 2 )(3 + 5p 2 )T l Tr 

AP (t)ap (3 + 10^ + 3^)^+2(1-^)^- ^ > 

Here, we used the same notation k as Ref.— to ease the comparison. Observable quantities, i.e., the charge, current, 
or spin averaged over measurement time, are either a symmetric (even) or an antisymmetric (odd) function of p 
depending on their symmetry properties under spin reversal. One can verify, for example, that J z averaged over the 
measurement time is an odd function of p. 

J zP = 0, 

j _ {J z )ap 2 P {4(i+p 2 )r L + (i-p 2 )r fl } 

zAP (t)ap (3 + 10p 2 + 3p 4 )r L + 2(l-p4)rV ^ °> 

Note that the average J z vanishes for P-alignment. 

At second order, one finds expressions, such as (/J2 is the centered moment of second order) 

MQ\p 



(i- P 2 )(3r L + 2r fl )3' 

8(1 - P 2 )(3 + 32p 2 + 38 P 4 + 40p 6 + 15p 8 )r £ r fl 
MaP - {(3 + KV + 3^ + 2(1-^)3 ■ (B6) 

Note the these correlation functions are actually cumulants or fluctuations around the mean value. At this order, 
expressions start to be lengthy, so that we list here only a few examples of our results: 

rn r L r fl {27(i + 3 P 2 )r| + 48 P 2 r L r fl + 4(3 + P 2 )r| 



(i- P 2 )(3r L + 2r K )3 

(1~P 2 )T l Tr(3T l + 2T r ) 



2(er 2 +4r L r fl + r 2 fl ) 

fM-JziAP - -r, 2vF -r tiT — t^tft ( B ^> 
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In order to check the consistency of these results, let us compare them with the FCS generating function Q(£, 77). 
For non-magnetic electrodes (p — 0), the analytic expression for £!(£,7y) is given in Refpi asji£ 



Z-zr L -r R , 1 



(B9) 



where z = (25 + 2)/(25 + 1), i.e., z = 3/2 for 5 = 1/2. Taking derivatives of Eq. (|B9[) with respect to counting fields 
£ or 77, one can, in principle, obtain any correlation function associated with Q and /, i.e., 



Km,n[Q,I}= lf((Q m n)c 



p=0 



1 d m d n 



(BIO) 



At lowest orders, this gives, 

ki[Q] p =o = Qp=o — 

K 2[Q]p=0 = V>2{Q]p=0 



3r L 



Kl[I]p=0 — lp=0 



3r L r R 

3Il + 2T R 



(3r L + 2r K ) : 



k 2 [/] p= o — /-*2[^] P =o — 



3T L T R (9T| + 4T|) 

(3r L + 2r fl ) 3 



(Bll) 



One can, therefore, check the consistency between Eqs. (IB3IB4IB6IB7I) and Eq. (|B9[) . by verifying the formulas given 
in Eqs. (|B"TT|) . The formulas (|B3IB4IB5IB6IB7IB8|) are also listed in Table I in the limit of T L -> 1 and T R — > 1. 
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TABLE II: The analytic formulas for the case of a molecular spin 5=1, when the state J = 3/2 of the occupied dot lies in the 
bias window. The first column shows the type of polarizations in the leads : P for parallel (Pl = Pr = p), AP for anti-parallel 
(Pl — —Pr = p), LP (Pl — p, Pr — 0) for left lead polarized only and RP (Pr = p, Pl = 0) for right lead polarized only. The 
second column shows the quantity whose analytical formula is given in the third column. 



